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Abstract 

Recent advances in modeling and control of crowds of pedestrians are briefly surveyed in this paper. Possibilities 
of applying fractional calculus in the modeling of crowd of pedestrians have been shortly reviewed and discussed from 
different aspects such as descriptions of motion, interactions of long range and effects of memory. Control of the 
crowd of pedestrians have also been formulated using the framework of Cyber-Physical Systems and been realized 
using networked Segways with onboard emergency response personnels to regulate the velocity and flux of the crowd. 
Platform for verification of the theoretical results are also provided in this paper. 

Keywords: Model for crowd of pedestrians, Control of crowd of pedestrians, Cyber-Physical Systems, Fractional order 
Calculus. 


1 Introduction 

As the most socially complex animals on the planet, research related to crowd of pedestrians has received a lot of attention 
in recent years. A lot of work has been conducted from perspectives of behavior, psychology, cognition and network to 
analyze problems related to particles, vehicles, robots, animal and even human beings (See [T], [2], [3], 0, [5]). Some 
manuscripts have been published in recent years concerning the modeling and control problem of crowd of pedestrians 

such as i, 0,0,0, m, eh, eh, eh, m, eh- 

On the other side, catastrophic events occurring around the world have demonstrated the need to re-analyze and 
re-examine current evacuation policies and procedures for crowd of pedestrians. The dynamic and uncertain nature of 
disasters has required the need of changing backup contingency plans according to evacuation needs. The first problem 
that confronted in the research of crowd of pedestrians is how to obtain a satisfactory model to characterize the complex 
nature of this kind of dynamic systems. There are a lot of characteristics which should be considered in the modeling and 
control problems such as self-organization, following leaders, common motives for action, psychological unity, emotional 
intensity, level of violence as shown in m- Due the features mentioned above, modeling and control of crowd of pedestrians 
are challenging tasks as shown in recent work of ei, ng, Ei, m, m and 1221 . 

Recent studies on modeling and control of crowd of pedestrians have been firstly reviewed in this paper in three 
different scales. Then necessities of introducing fractional calculus in the modeling and control of crowd of pedestrians 
have been discussed in detail and networked Segways with on-board emergency response personnels have been employed 
in the evacuation control of crowd of pedestrians. Framework based on Cyber-Physical systems have been formulated in 
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the end for studying the modeling and control problem of crowd of pedestrians. Main contributions of this paper lie in 
that the tool of fractional order calculus have been introduced in the modeling and control of crowd of pedestrians where 
some phenomena observed from realities are impossible or not easy to be considered just in the framework of integer order 
calculus. The aim of this paper is to provide an alternative way for studying the complex crowd of pedestrians which is 
much more close to reality. 

The rest of the paper is organized as follows. Model for crowd of pedestrians are firstly reviewed in three different 
scales in Section [2] respectively; Then recent advances in control of crowd of pedestrians are discussed in micro-scale and 
macro-scale in Section [3] Our framework for studying the evacuation problem of crowd of pedestrians are shown in Section 
Q] based upon calculus of fractional order. Validation of obtained theoretical results are also included in this Section; 
Conclusions of this paper are given in Section [5] 


2 Modeling of Crowd of Pedestrians 

A lot of modeling methods have been proposed recently. Most of these modeling methods are based on inspirations 
that the crowd of pedestrians can be treated as moving particles, smoothing fluids, or granular flows in different scales. 
Modeling methodologies are firstly reviewed in micro-scale, macro-scale and mesco-scale, respectively. Then necessity of 
introducing fractional order calculus in the modeling of crowd of pedestrians are discussed and some recent advances of 
modeling using fractional order calculus are also included in this section. 


2.1 Micro-scale—Ordinary Differential Equation (ODE) 


When the density of pedestrians is low, each pedestrian can move freely and interactions among pedestrians can be 
modeled using the framework of social forces. A lot of research has also been conducted using the Newton’s laws as shown 
in the following microscopic model 
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where ff is the self-driven force towards some desired velocity, /A is the interaction between agent % and its neighbor j 
and j'lf represents the interactions with environment such as walls or corridors. 

Different and ff v have been constructed in the work of [23] for panic scenarios and it has been proved that neither 
internal forces nor external forces alone can perform well in the evacuation of crowd of pedestrians and the practical way 
is to consider them simultaneously in the micro-scale model ©; Similarly, it has been shown in |24j that not only external 
forces such as interactions with neighbours and environments but also some internal forces such as will force and personal 
information should be considered in the modeling process. But the fact is that a lot of work has been done from the 
view of multi-agent system where only external forces are considered. Interactions with external environment such as the 
spatial-temporal markings have been described as potential function and added into the right-side of m as done in [25] ; 
Relationship between formation patterns, stability and types of interactions have also been analyzed in under external 
potential function to find whether a swarm will collapse or not. Some recent work such as m, m\ from the aspect of 
modeling and |29j . [30] from the aspect of control of crowd of pedestrians are also only focused on external interactions 
with neighbouring agents. 

Most of the simulation results have been conducted in this scale as shown in [12] because the methods are simple and 
animations are realistic. In order to satisfy the requirements of high realism and real-time animation in the modeling 
and control in this scale, the idea of mapping desired behaviors to the stable solutions of classical nonlinear dynamic 
system such as Van der Pol oscillator or fix-point attractor has been introduced in [31j and [32l where one nonlinear 
transformation connecting pedestrian’s periodic or non-periodic motion with structurally stable system has been firstly 
constructed and then real-time animation problem can be solved in low-dimensional space without losing details; Methods 
from cognitive science such as heuristic behavior have been employed to adjust the walking direction and speed in micro¬ 
scale in pQ; Global roadmap-based navigation method has been proposed in the homing or exploring problem as shown 
in [33] . Although these methods are simple and efficient, the prescription of all behaviors is not very easy for crowd of 
pedestrians of large numbers. 

Remark 1. Comments on Microscopic Model: 


• Main advantages of the microscopic model lie in that heterogeneities of each pedestrian can be considered explicitly 
and simulation results obtained are of highly realistic. But the microscopic model is not a good choice if the number 
of pedestrians are very high as some unnecessary interactions or effects have been included; 

• Most of the previous work has just treated each pedestrian as one physical particle and few works have been done 
under the consideration of pedestrian’s memory or some other internal effects; Local interacting rules have received 
a lot of attentions in previous research while long range interacting rules are not so popular. 
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2.2 Macro-scale—Partial Differential Equation (PDE) 

The motion of crowd of pedestrians show some striking analogies with the motion of fluids when the density goes high. 
Thus research of crowd of pedestrians in macro-scale has benefited a lot from the macro-scale research of traffic system 
where the well known LWR (Lighthill-Whitham and Richards) model d2j and PW (Payne-Whitham) model © have been 
proposed. It is well known that the traffic dynamic system in macro-scale as shown in m can be described by 

= ° ( 2 ) 
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where p denotes the density of crowd and f(t,x) is the flux of crowd, V{p) is equilibrium Speed, v is the relaxation 
term, A ^ x is the anticipation term and //A 221 is the viscosity term. Relationship between flux and density has attracted 
a lot of interests and a lot of work has been done to generalize the LWR model and PW model. 

Two versions of LWR model such as high-order LWR model and low-order LWR model have been shown in [Ml . 
Compared with classical LWR model that only concerns transitions among equilibrium states and homogeneity, both 
non-equilibrium transitions and inhomogeneities have been accommodated in the work of [34] through constructing much 
more complex relationship between flux and density; Similarly, new relationship between speed and density has been used 
to derive Lighthill-Whitham model in panic scenarios in [35]. 

Also based on the fundamental conservation law of mass and momentum, macroscopic model for crowd pedestrians in 
two dimensional space, as below 

|„+Am + Ato) = o (4) 

has been constructed in [361 for different types of pedestrians with different walking habits. Macroscopic model in two 
dimensional space has been generalized in EZ3 with terms for anticipation and relaxation. 

As interactions among pedestrians are assumed to be the same in macro-scale therefore it is not easy to incorporate 
the heterogeneity of each pedestrian. In order to characterize the crowd of pedestrians more precisely, models that can 
benefit from both the micro-scale and macro-scale are much preferred. A modeling procedure based on the time-varying 
measures 

pt = 0m t + (1 - 0)M t , (5) 

has been prescribed for the multi-scale characterization of crowd-pedestrian system where m* = 4-P ? (t) and dM t {x) = 

pit , x)dx are the microscopic and macroscopic mass respectively. Both topological interactions such as the Braess’s paradox 
(obstacles may contribute to improve the flow of people in some situations) characterized in macro scale and anisotropic 
interactions such as the granular role of some pedestrians described in micro-scale can also be considered in this framework 
through the time-varying measure © ; Similar results that are composed of agent-based microscopic model and flow- 
based macroscopic model have been shown in ©I and [321 where the accuracy of microscopic model and the efficiency 
of macroscopic model have been combined together through initializing each other by techniques of aggregation and 
disaggregation. 

Remark 2. Comments on Macroscopic Model: 


• Computation time have been greatly decreased as each pedestrian has been treated as physical particle and have 
same characteristics in macro-scale; 

• Main disadvantages of macroscopic models lie in that heterogeneities of pedestrians such as interactions or mobilities 
cannot be characterized or considered in detail in this scale. Obtained results using the macroscopic model just 
provide some references for the control of crowd of pedestrians. 


2.3 Mesco-scale—Integral Differential Equation (IDE) 

After Boltzmann formalized the concepts of kinetic equations in the nineteenth century, it has been widely used in 
astrophysics, engineering, social science and even in biology. Although the following basic kinetic-transport equation 

+£ • V x /(t,ar,^) = 0 (6) 
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that describes the evolution of density f(t, x, f) of particles passing through (t, x ) with velocity f seems very simple, great 
success has been achieved in the applications of Boltzmann’s kinetic equation and many generalizations of this equation 
have been made. 

A lot of modeling results for crowd of pedestrians have been derived from the equation © since macro-scale variables 
and micro-scale variables have combined with each other in this equation. Nonlinear integral-differential equations have 
also been used in [3D] to describe the crowd of pedestrians with competitions and kinetic interactions as 

df 

= Ji[f}(t,u) + 7 , (t, u) 

where the Ji[f](t,u) describes the evolution of density due to gain and loss of distribution function and 7 i(t,u) describes 
the production and migration of the group; Mescoscopic models with binary interactions, averaged binary interactions 
and mean field interactions have been given in [41] as 


df df d(fF[f]) 
dt ^dx df, 


Q(f,u) 


(7) 


where different F[f] and Q{f,u) are employed to describe these interactions. Conservation law of mass and momentum 
in macro-scale has been generalized in [3^1 where internal intentions and external interactions have been simultaneously 
considered in the following mescoscopic model 


d_ d_ 

dt Pti + dx 
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where p^ is the density, the first term on the right is caused by pedestrians entering or leaving some interested areas 
and the second term describes the effects caused by internal intentions and external interactions. Equation (JHJ has been 
generalized to the following two dimensional mescoscopic equation 


dtp + [d Xl (pv 1 ) + d X2 {pv 2 )] + [d Vl {pAi) + d V2 {pA 2 )\ 

= (dtp) event ~ {dtp) event 

in |43] where d Xl {pv 1 ) + d X2 {pv 2 ) means the changes of density due to convection, d Vl {pA{) + d V2 {pA 2 ) are the terms of 
acceleration and deceleration, {dtp)+ vent — {dtp)f ven t means the interaction of events. 

Some other research have also been conducted from the aspect of micro-scale using the framework of mesco kinetic 
theory. With the help of activity variables, heterogeneities of each pedestrian have been characterized in micro-scale in 
[44] and [45] using the statistical distribution of position and velocity; Enskog-like interactions and stochastic interactions 
have been considered in [46] ; Nonlinear interactions instead of linear interactions have been applied in [47] ; Both short and 
long range interactions have been reported in [44] contrary to previous modeling methods where only local interactions 
are assumed. 

Remark 3. Comments on mescoscopic model: 

• As both information from micro-scale and information from macro-scale have been included in this framework, 
heterogeneities of pedestrians and different kinds of interactions can be included easily in the right hand side of ©; 
Another advantage is that efficiency of the mescoscopic model has been greatly improved compared with that of 
microscopic model; 

• One of the main disadvantages of mescoscopic model is that existence of analytic expressions for equilibrium is not 
guaranteed for obtained IDEs now and another disadvantages of this framework is the technical difficulty when there 
are more than one microscopic variables that needed to be included in the above IDEs; 

2.4 Fractional Model 

Fractional Calculus has shown great potential in different applications such as particles in fluid, plasma physics, quantum 
optics and many others. Some phenomena such as self-similarity, non-stationary and spiky phenomena, non-Gaussian 
relaxation, short or long memory and long range interactions are all closely related to fractional calculus. The authors have 
noticed that previous research has mentioned or implied potential applications of fractional calculus in modeling, control 
and optimization of crowd of pedestrians in different aspects such as description of motion, description of interactions and 
hysteresis phenomenona. Recent advances in these aspects are shortly reviewed in the followings. 
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1. Fractional Calculus in Description of Motion 

Some applications of fractional calculus has been done for physical particles in macro scale. Based on the fractional 
derivative of order a in [481 


d a f(x,t) = f(x,t + At) - f(x,t) 
dt a A™o A t a ’ 

some different fractional kinetic equations (FKE) and their properties have been explained in [5S] such as the FKE 
with desired direction, anisotropic FKE and nonlinear FKE; Then fractional kinetics of distributed-order has been 
proposed in [55] and it has been shown that the fractional model with distributed order can be effectively used to 
describe the accelerating and decelerating of sub-diffusion and super-diffusion. 

Based on these preliminary works, fractional model for crowd of pedestrians in the macro-scale has been proposed 
in m using the fractal derivatives of time and space where fractional order dynamics for one dimensional crowd 
can be described as 


d d 

x ) + X )1 = 0 ( 10 ) 

where p(t, x ) and v(t, x) are the density and velocity of crowd at position x and time t, a and /? is the fractal order 
of time and space respectively. Similar results can also be found in [52] where a = /3 has been assumed. Considering 
that each pedestrian’s decision making process plays an important role in his (her) next movement, the dynamic 
process of decision making has been described in [21] where coupling dynamic models composed of the Fractional 
Fokker-Planck equation 


d a P(u, t) _ d[f(u, d, t)P(u, t)\ , ^ 

dt a - du +~{PU,t) [Li) 

and the Fractional Master equation 

d = d, t) + g{u , d, i)£(i) + rj(t) (12) 

have been constructed and P(u,t) is the distribution function, E(P, u,t) is the additional terms due to variations 
in infinitesimal time increment, u(x,y,t ) defines the cost that each pedestrian should pay for desired destinations, 
f{u, d, t) is the nonlinear dependence on previous costs u(x, y, t ) and previous decisions d(t). The forward dynamics 
m describes the evolution of distribution function with respect to time while the backward equation m describes 
the evolution of pedestrian’s decision with respect to the inverse of time. In order to control and predict the evolution 
of crowd of pedestrians, optimal control of this complex system has been firstly reformulated in [53] based on the 
fractional master equations. 

Concerning the fractional model in micro-scale, fractional Langevin equations have been given in [54] for human 
motion tracking and fractional models with fixed order and variable order have been proposed in [55] and [56] to 
describe motion of pedestrians. 

2. Long Range Interactions of Power Law 

Different formation patterns for swarming of fish can be generated through tuning the area of attraction and repulsion 
as shown in El; Quantitatively characterizing of relationship between the area of repulsion and attraction have been 
analyzed in [58] for dynamic model of integer order. Thus even for dynamic models of integer order, phenomena of 
transiting from ordered state to disordered state have been observed in [55] under long range interactions. Transition 
from chaos to turbulence has been reported in m where fractional long range interactions of (l is the distance 

between oscillators, a is the tuning parameter) has been imposed for nonlinear oscillators. Recently, transitions 
have also been declared for Boltzmann equations using local, non-local or configuration-dependent interactions in 
m where transition from purely diffusive regime to flocking patterns can be realized by modulating the range and 
strength of interactions. We can say that range of interaction and power of interaction play an important role in 
study of system of large numbers. 

Actually, the long range interactions are connected with fractional dynamics. Relationships between long range 
interactions of power law in micro-scale and fractional Euler-Lagrange equations in macro-scale have been discussed 
in [52] and it has been proved that there exists one fractional Euler-Lagrange equation in macro-scale if there are 
long range interactions of power law in micro-scale; From which we can say that the gap between micro-scale and 
macro-scale can be bridged together using long range interactions of power laws. As the long range interactions are 
so common in crowd of pedestrians, it is reasonable to describe the dynamics of crowd of pedestrians using fractional 
calculus. 
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3. Hysteresis phenomenon and Collective Memory 

It has been shown in E3 that minor changes in individual’s responses can lead to different collective behaviors such 
as alignment, swarm and torus where the collective memory has played an important role; Hysteresis phenomenon 
has been observed in [1] where nonlinear relationship between collective behaviors and range of interactions has been 
explicitly shown in Figure 1 of {33. 

According to the following Definition of Hysteresis from Wikipedia: Hysteresis is the dependence of the output of 
a system not only on its current input, but also on its history of past inputs. In other words, the change of 
group’s behavior not only depends on the current control input but also depends on the history of individual behavior 
and the shape of group; The hysteresis phenomenons or effects of short memory are so common that they can be 
found in the acceleration, deceleration and equilibrium of traffic flows in |63| and individual’s actions such as buying 
and selling in the stock market. For particles with long time memory, fractional Fokker-Planck equations have been 
derived in [64] based on the correlating functions of probability densities that follow power law; A non-negative 
memory function ?](■) has been introduced in [63] to generalize the LWR model using the form of convolution 

P*V= r}{x - £)p{t,£)d£ 

Jr 2 

where preferred path can be found and regions of high density can be avoided. 

As calculus of integer order can be considered as special case of fractional order calculus, each pedestrian can be 
characterized much closer to reality using the framework of fractional order calculus. Many effects such as individual’s 
memory, long-range interactions which are difficult to characterize using calculus of integer order can be compensated 
using calculus of fractional order; 

Remark 4. Compared with dynamic models of integer order in previous research, some advantages and disadvantages of 
dynamic model of fractional order are listed in the followings: 

• In time domain: Only normal diffusive process has been considered in previous study in macro scale due to the 
limitation of calculus of integer order. Besides the normal diffusive process, sub-diffusive process and super-diffusive 
process have been added to describe the crowd in macro scale using fractional calculus in time domain; 

• In spatial domain: Dimension of space are only limited to 1, 2 or 3 in previous research while the fractional 
dimension can be included in fractional calculus to consider effects of environments on crowd pedestrian system; 

• Some other elements: As we have shown that some phenomenons such as self-similarity, non-stationary and spiky 
phenomenons, non-Gaussian relaxation, effects of memory, long-range interactions which are difficult to characterize 
or explain using calculus of integer order can be compensated using calculus of fractional order; 

The authors do not want to prove that previous models of integer order are not effective anymore in reality and we 
just want show that Fractional Calculus as generalization of Calculus of Integer order has provided us much more freedom 
in characterizing and understanding of the complexities of crowd of pedestrians. The authors admit that there are a lot 
of challenging work left to do for the obtained model of fractional order such as controller design, stability analysis and 
performance evaluation. Most of them are not so easy to solve now. 


3 Control of Crowd of Pedestrians 

3.1 Control of Macroscopic Model 

In the evacuation of crowd-pedestrian system, moving direction and moving speed are of great importance to guarantee 
the smooth movement of crowds. Some previous research has chosen v(t,x) in ([3]) as the control input and distributed 
state feedback controllers have been constructed in [66] and [671 for the macro-scale model in one and two dimensional 
space. Main idea of [BB| and m is using finite dimensional system to approximate infinite dimensional system; Then 
theory of nonlinear control can be used for design of feedback controllers for this approximating system. Main problems 
of this framework is that the original system may be still unstable even if the obtained controllers work very well on 
the approximating system. Different from the framework of approximation, the control and stability problem of the 
macroscopic model ([3]) has been formulated directly in the distributed control of partial differential equations in [68] and 
[69] . Previous controllers have been generalized to diffusion, advective and advective-diffusion controllers. Comparisons of 
these controllers have been done in EZQ! and it has been shown that much more serious problems may arise if only diffusion 
feedback controller is adopted as there is no preferred direction for the pedestrians to follow in evacuation process; But 
faster evacuation can be realized for LWR model under diffusion-advection-state feedback controllers since the direction 
of evacuation is provided. 


6 


In order to avoid undesirable congestion and blockages in the evacuation of pedestrians, optimal feedback controllers 
have been proposed in m and m by instructing pedestrians to adjust their velocities. Although optimal results have 
been obtained in m and m, applications of these results are only limited to the deterministic case. Robust control 
of crowd-pedestrian system has interested a lot of people in recent years. For example, Lyapunov techniques have been 
firstly utilized in m to construct velocity controllers for automated highway systems where not only position and time 
but also effects of lanes, drivers and destinations have been included in the macroscopic PDE model; Previous diffusion- 
state feedback controllers have been generalized using methods of Lyapunov redesign in [73] to deal with disturbances in 
control input; Robust controllers based on sliding model control have been constructed in m for the control of crowd of 
pedestrians with matched and unmatched uncertainties that are caused by external disturbance and parametric variations. 

Remark 5. Similar to the modeling methods for macroscopic model, heterogeneities of controllers have also been neglected 
in this scale. 

3.2 Control of Microscopic Model 

3.2.1 Control using Leaders with Naive followers 

Based on inspirations from swarming of ants, schooling of fish and flocking of birds, how to formulate or control the 
collective behaviors in micro-scale has received a lot of attentions from the community of control, computation and 
computer. One of the fascinating phenomenon is that this kind of collective behaviors can be generated from simple or 
basic interacting rules and control of this kind of collective behaviors can be realized through controlling just a small part 
of the group. 

Determining of moving direction based on neighbor’s direction have been firstly proposed in m for average consensus 
problem of each agent. Work of [5] and [3] has explicitly shown the quantitative relationship between individual’s local 
interacting rules and group’s collective behaviors; Three experiments have been done in m to show the effects of number 
and topology of informed agents on the consensus problem and it has been proved using experiments that only a small 
number of informed individuals are enough for driving a large number of uninformed individuals to reach consensus 
without explicit communications; Further results concerning consensus under uncertain or conflicting information have 
been considered in [78]. For evacuation problem of crowd pedestrian system, result of tzg has shown that the evacuation 
rate can be greatly improved by adding some leaders with global knowledge and it is always good for the evacuation process 
if we place some leaders in the immediate proximate of crowds and some leaders scattered around the environment. 

Remark 6. Similar to the pinning control of complex network with large number of nodes, how many nodes are needed to 
control the network and which one should be chosen as the pinning node are interesting problems. For crowds, how many 
leaders and what kind of leaders are needed in control of crowd of pedestrians are worthy of further considerations. 

Phenomena such as self-similar structure, heavy-tailed workloads and long-range dependence effects that have been 
proposed in m, es for study of network are related to calculus of fractional order. The authors think that not only 
topology and connectivity should be considered but also the dynamic evolutions of the network itself should be considered 
especially when confronted with stochastic network of large numbers Thus control in the framework of fractional order 
calculus is very interesting although there are many challenges. 

3.2.2 Control without leaders 

• Decentralized Framework: 

Using the methodology of decomposition, control problems for complex system can be transformed into control 
problems for simpler subsystems. A lot of decentralized controllers have been obtained under the framework of 
decomposition. Decentralized controllers for large number of stochastic agents have been given in (821 where not 
only the evolution in time scale but also the evolvement in “space” scale ( where TV —>- oo, TV is the number of 
agents) are considered; Decentralized controllers have also been constructed using the framework of decomposition 
in (83] where complex LQG games problem has been reduced to two-player games problem. 

Advantages of the decentralized framework are that the burden of computation has been greatly reduced and obtained 
results can be easily extended to systems with large numbers. However, neighbor’s information is not used in the 
decentralized framework, where it is much preferred to have distributed controllers. 

• Distributed Framework: 

Some other distributed feedback controllers for macroscopic model have been reviewed in previous sections or can 
be found in m, E3 and [f)S] ■ Some consensus protocols of multi-agent system also belong to this framework which 
we will not state one by one and interested readers can find some recent advances in the review papers of [811 . [85] 
or references in them. 

One big challenge for distributed control lies in the burden of computation if the number of pedestrians goes to 
infinity. Some previous research has adopted the mean field methods to estimate the influence of neighboring agents 
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to relieve the burden of computation. Mean field methods have been used in [BB; for controller design for leader- 
follower agents with linear stochastic dynamics; Distributed controllers for agents with one major/leader agent have 
been given in m with Markov jump parameters in controlled system and random parameters in objective functions 
where mean field method has been utilized to estimate effects from neighboring pedestrians and the leader; Mean 
field LQG controllers have been proposed in [88] for socially optimal control problem whose cost functions are coupled 
with each other; Mean field control strategies have been given in [85] for the consensus problem of multi-stochastic 
agents based on the assumption that information of initial state distribution is available to each agent; Similar results 
can also be found in [90] in the control of stochastic multi-agent systems where statistical information of neighboring 
agents are obtained using mean field methods. 

In order to study the true interacting rules among agents, game theory has also been combined with mean field 
method in control of agents of large numbers which we will call it mean held games in the followings. Based on 
mean Held games introduced in [9T| . dynamics of human’s decision making process has been considered in [921 
where mean Held game theory has been employed in the backward Hamilton-Jacobi-Bcllman (HJB) equation and 
forward Fokker-Plank (FP) equation; The Forward-backward equations have also been adopted in [93] to describe 
the aversion and congestion phenomena in macroscopic scale where each pedestrian’s ability of anticipation has been 
characterized using the backward equation; 

Advantages of the mean field method is that only information of initial distributions is needed in the design of controllers 
and the popular topology condition such as connected graph or jointly-connected graph is no longer needed; Additional 
benefit is that communications among neighboring individuals are no longer needed and obtained results can also be easily 
extended to systems with large numbers. 

Remark 7. Research of crowd of pedestrians has benefit a lot from the research of traffic control system where limiting 
speed has served as an effective way to guarantee the normal flow of vehicles. Information of density has been used in 
specifying limitation of speeds in different area in |14j . Roadmap-based planning has been used in [28] for evacuation of 
large number of agents where interactions and coordination among agents are of great importance in realization of the 
evacuation. 

Speed, density and even interactions are the main information that has received a lot of attention from the view of 
multi agent system. Besides the above information, it is the authors’ belief that effects of memory, environment or even 
building structure should be considered in modeling and control of crowd of pedestrians where Fractional Calculus will 
play an important rule. 


4 Our Framework 

A fractional framework for modeling and controlling of crowd of pedestrians has been proposed in this paper. The proposed 
research represents the coupling of the ‘PHYSICAL PART’ (Mass Pedestrian Evacuation management of crowds) with 
the ‘CYBER PART’ (modeling and prediction of crowds). And transferring of information between these two parts has 
been implemented through networked Segways, with on-board emergency response personnels, and facility sensing and 
actuation. 


CYBER PART 


Modeling of Evacuation 


Ordinary Differential Equations 


Partial Differential Equations 
Integral Differential Equations 


Coupling Equations 



PHYSICAL PART 


Controlling of Evacuation 


Movement of each Pedestrian 


Smoothing Fluids of Crowds 
Granular Fluids of Crowds 


Actions Based on Games 


Figure 1: A Cyber Physical Framework for describing the coupling of model, control, and management 

In Cyber Part, ordinary differential equations, partial differential equations and integral differential equations have 
been employed to describe the crowd of pedestrians using calculus of fractional order in micro-scale, macro-scale and 









































meso-scale, respectively. Interesting information such as speed, density, flux and even formation patterns that obtained 
through CCTV, Segways, Cellphone and some other sensors can be used in calculating the obtained models, controlling 
the crowds and even predicting the stampede that is going to occur. Closed loop system that are composed by the cyber 
part and physical part in Figure [l] has shown the flow of information and major points of each part. 

4.1 Fractional Modeling of Crowd of Pedestrians 

• In Micro-scale level with low density of crowds, behavior of each pedestrian can be described by the ordinary 
differential equations based on some widely used methods such as social force model © or agent based model where 
parameters of systems can be further calibrated using empirical or observed data. Main reason of using these classical 
models is due to the convenience of introducing heterogeneities in the description of pedestrians in micro-scale; 

• In Macro-scale level, the density of crowds is so high that the motion of all pedestrians can be modeled as continuum 
fluids where partial differential equations can be derived from the conservation law of mass or momentum on 
interested area. Main differences between crowd of pedestrians and smoothing fluids are that different motion 
patterns such as crossing or intersecting is allowed in crowd-pedestrian system due to the freedom of choosing 
different routes. The authors believe that both the time scale and the spatial scale should be considered in the 
modeling of crowd of pedestrians and preliminary work of fractional macroscopic model (1101) has been shown in [51] 
where the fractional order associated with time and the fractional order associated with the fractal structure have 
been included; 

• In Mesco-scale level with medium density, the dynamics of evacuation or egress process is similar to the diffusion 
process of active particles in many aspects such as the porous or granular pattern in the smoothing fluids. Fractional 
convection-diffusion equations are useful tools to model the crowd pedestrians in this case due to the phenomenon 
of porosity observed in this scale; Heterogeneous pedestrians can be modeled using mobile potential fields indexed 
by activity variables as shown in Figure [2] which can be included in the convection-diffusion equations to guarantee 
the existence of heterogeneities on this level; Interactions between microscopic model and macroscopic model can 
also be realized through mean field games to increase the validation of obtained models and relieve the burden of 
computation. 



Figure 2: Normalized effect of activity variables on potential fields (PFs) 


Reasons of using the above modeling framework: 

• As there is no general method prescribed for modeling crowd pedestrians for all kinds of scenarios, it is much 
reasonable to choose the most appropriate model for different problems. Since the microscopic model is powerful to 
describe the heterogeneities of pedestrians, we choose to use the microscopic model when the density is low; With 
the increasing of density, granular flows with porosity phenomenon can be observed in large crowds. This kind of 
heterogeneities are modeled using different mobile potential fields which can be easily added into the right hand side 
of © to describe the effects from micro-scale; As density grows, the porosity or granular phenomenon of smoothing 
fluids disappears and fractional dynamic model can also be obtained by using generalized conservation law of mass 
or momentum as deduced in previous research. 
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• The macroscopic model is responsible for generating a homogenizing effect with desirable smoothness in macro-scale, 
and the microscopic model is responsible for characterizing heterogeneities and interactions from the macro-scale. 
For mescoscopic model, not only heterogeneities and porous patterns of crowds can be explicitly characterized but 
also interactions between micro-scale and macro-scale can be included in the obtained integral differential equations. 

• Although we describe the modeling of crowds in different scales according to their densities, these obtained models 
are not independent of each other as shown in Figure[3] The macro-scale variables such as density or flow comes from 
aggregation of micro-scale data of each pedestrian and the motion of each pedestrian is also affected or constrained 
by the people around him. 



Figure 3: Fractional Model in Micro-scale and Macro-scale 


4.2 Fractional Control of Crowd of Pedestrians 

• In microscopic level, the control of each pedestrian is focused on analyzing the relationship between different interac¬ 
tions and collective patterns where not only interactions based on topology but also interactions based on short range 
and long range interactions will be considered. Dynamics of crowd-pedestrian system under long range interactions 
are closely related to fractional calculus and have an important effect on generation of different collective patterns; 

— This phenomenon observed from self-organized groups under different interactions can be used to control the 
collective motion of crowds where mobile Segways have been added into the group to tune the range interactions 
among pedestrians to generate some desired collective patterns; 

• In macroscopic level, fractional controllers based on fractional convection, diffusion or both of them will be con¬ 
structed for controlling of the smoothing fluids m for crowd of pedestrians with fractional time/spatial orders; 

— Due to the high densities in this scale, it is not easy to inject control agents into the fluids. In our framework, 
Mobile Segways with emergency personnels will be dispatched to control the inflow and outflow of crowds from 
outside based on the theory of boundary control to guarantee the smoothing evacuation of high-density crowds 
without breakdowns; 

• In mescoscopic level, forward fractional convection-diffusion equations and backward fractional Hamilton-Jacobi- 
Bcllman (H-J-B) equations are used for modeling the crowd pedestrian dynamics where the forward part describes 
the evolution of crowds and the backward equation describes the evolution of decision making process. 

— Movement of next step will be prescribed based on the probability distribution function as shown in backward 
part (fl2ll of the fractional HJB equations; 

— Fractional order controllers based on fractional diffusions or fractional convections are adopted to realize the 
control or optimization of the forward part ( 1111) of the fractional HJB equations; 

— Mean field games will play an important role in estimating neighbor’s influence and reducing burden of com¬ 
munication and computation; 


10 































— Mobile Segways with global instructions can be used to guide or drive the crowds through broadcasting or 
changing the structure of space to control the velocity or flow of the crowds; 

Reasons of using the above control framework: 

• For crowds with low density, mobile Segways can be easily added into the crowds and control of the crowd-pedestrian 
system can be realized through interactions with just a few of them to realize control of each pedestrian; 

• For crowds of high density, it is not wise to add mobile Segways into the crowds again due to the high density. In 
this case, movements of pedestrians are firmly constrained by the people around them. What we try do is using 
more delicate controllers such as boundary control to tune the inflow and outflow of people through Segways to 
realize the smoothing fluids without breakdowns; 

• For crowds with medium density, there are some concerns that are not easy to be solved in this scale. Besides 
heterogeneities and computation burdens, the interactions between microscopic model and macroscopic model are 
the main challenges. In other words, balancing between realization of each pedestrian’s target and desired evolution 
of the entire group can be manipulated simultaneously in the framework of dynamic game theory based on mean 
field. 

4.3 Platform for Verification 

Since real verifications of the theoretical results for evacuation of crowd of pedestrians are not easy to conduct, a lot of 
simulation results have been done in previous research. There are a lot of software for the simulation research such as 
VISSIM, EXODUS, Simulex, PSCrowd, PEDSIM, and VISWALK, et al. 

• VISSIM is one powerful tool available for simulating multi-modal traffic flows, including cars, buses, motorcycles, 
bicycles and pedestrians and it is also a useful tool for the evaluation of various alternatives based on transportation 
engineering. In VISSIM, movement of each pedestrian is modeled by the social force model ([94]) where a total 
force resulting from the social, psychological, and physical forces has been imposed. The forces that are influencing 
pedestrian’s motion are caused by his/her intentions to reach his destination as well as by other pedestrians and 
obstacles. Thereby the other pedestrians can have both attractive and a repulsive influences. 

• Due to the requirements of data processing and computations in modeling and control of crowd of pedestrians, 
the simulation research is much preferred to be realized using Matlab. The authors of this paper try to study the 
modeling and control problem using the platform of DIFF-MAS that has been developed in [55] for simulating the 
measurements and control of diffusion processes using Matlab script and Simulink. Optimal placements of sensors 
and actuators, and structures for understanding distributed networked actuation and sensing have been done on this 
platform. Recently, this platform has been updated for measuring and control of diffusive process of fractional order 
which is useful for modeling and control crowds of pedestrians using calculus of fractional order. 

Remark 8. Considering unexpected or dangerous events in real-life experiment, only initial theoretical results and some 
simple simulation results have been done by the authors under the framework of fractional order modeling and control of 
crowd of pedestrians. There are lots of work left to show the effectiveness of this framework in future especially from the 
aspect of control. 

5 Conclusion 

Recent advances in modeling and control of crowd of pedestrians have been reviewed in the framework of Cyber-Physical 
Systems. In the first part, not only modeling methods in micro-scale, macro-scale and mesco-scale are reviewed but also 
possibilities of applying fractional calculus are discussed; Due to challenges existing in the control problem, only controllers 
in macro-scale and micro-scale have been surveyed such as state feedback controllers, distributed controllers and robust 
controllers. Initial considerations about control of crowd of pedestrians are discussed in the framework of Cyber-Physical 
Systems. Long range interactions are used to generate self-organized collective behavior in the micro-scale; Fractional 
diffusion-convection controllers will be constructed for crowd fluids with fractal time-spatial orders in the macro-scale; 
Fractional Mean Field Games theory will be adopted in modeling pedestrian’s decision making process and control the 
fluids of crowd with porosity in meso-scale. Multiple mobile Segways with onboard emergency response personnels are 
employed to realize the control of velocity and flux of crowd of pedestrians in different scenarios. 

A lot of problems are still open and needing cooperation of multi-disciplines such as environmental design, engineering, 
and transportation. Potential topics such as modeling the crowd of pedestrians with consideration of stochastic noise, 
psychological effects and constructing controllers under requirements of scalability and robustness are both inter¬ 
esting topics and worthy of much efforts in future. Another interesting and important topic is concerning the security 
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problem in the crowd of pedestrians where some recent results on security of consensus problem such as m, m and [HE] 
are beneficial to solve of this problem. 
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